Azimuth determination for vector sensor tools

ABSTRACT

This invention relates to mapping or survey apparatus and methods, and more particularly concerns derivation of the azimuth output indications for such apparatus in a borehole from the outputs or output indications of either an inertial angular rate vector sensor (or sensors) and an acceleration vector sensor (or sensors), or a magnetic field vector sensor (or sensors), and from the outputs of an acceleration vector sensor (or sensors). At least one of such sensors in any instrument may be canted relative to the borehole axis. Borehole tilt is also derived.

BACKGROUND OF THE INVENTION

This application is a continuation-in-part of our prior application Ser. No. 351,744, filed Feb. 24, 1982 now U.S. Pat. No. 4,433,491.

This invention relates generally to mapping or survey apparatus and methods, and more particularly concerns derivation of the azimuth output indications for such apparatus from the outputs or output indications of either an inertial angular rate vector sensor (or sensors) and an acceleration vector sensor (or sensors), or a magnetic field vector sensor (or sensors), and from the outputs of an acceleration vector sensor (or sensors).

This invention also relates to methods used to compute tilt and azimuth from outputs of canted sensors which have a component or components of their input axis or axes of sensitivity along the axis of rotation, for survey apparatus having an axis or axes of rotation for one or more of the sensors.

U.S. Pat. No. 3,753,296 describes the use of a single inertial angular rate sensor, or "rate-of-turn gyroscope", and a single acceleration sensor, both having their input axes of sensitivity nominally normal to the direction of travel in a borehole and parallel to each other for survey in a well or borehole. In this case, both sensors are rotated about an axis parallel to the borehole by either the carrying structure and container or by a rotatable frame internal to the survey tool. U.S. Pat. No. 4,199,869 describes the use of one or two dual axis inertial angular rate sensors in combination with a dual axis acceleration sensor for survey in a well or borehole. Again in this case, the sensors are rotated about an axis parallel to the borehole by either the carrying structure or by a rotatable frame internal to the survey tool. U.S. Pat. No. 4,244,116 describes the of a one dual axis inertial angular rate sensor having its spin axis parallel to the borehole axis and one dual axis accelerometer for survey in a well of borehole. In this case no provision is made for rotation of the sensors about the borehole axis. U.S. patent application Ser. No. 338,261, filed Jan. 11, 1982 describes the use of one or more magnetic field vector sensors in combination with one or more acceleration sensors for survey with respect to the earth's magnetic field vector in a way related to the sensors of U.S. Pat. Nos. 3,753,296 and 4,199,869 which survey with respect to the earth's inertial angular rate vector.

The referenced patents and application describe the sensing equipments and show provisions to compute the output desired azimuth indication, but none of them show or teach the method and means described herein for obtaining the desired output, nor do they show the essential use of the output of the acceleration sensor (or sensors) to resolve the output (or outputs) of either the inertial angular rate sensor (or sensors) or the magnetic field sensor (or sensors) into a known coordinate system. For example, U.S. Pat. No. 4,244,116 shows a computation of: ##EQU1## Where E=Vector of the earth rotation

Ω_(z) =Component of E along the borehole

Ω_(x), Ω_(y) =Gyro outputs normal to the borehole,

and states "the measurement Ω_(x) and Ω_(y) and the calculation of Ω_(z) give then the azimuth of the drilling line". This is in general not true since Ω_(x) and Ω_(y) are known only to be perpendicular to the drilling line, but are not known in a known earth fixed coordinate set.

It may be shown that the apparatus of the above cited patents and applications, and with respect to described methods of computation, have associated geometric regions in which poor accuracy can result. Thus, for the determination of tilt or inclination, poor accuracy results when the plane containing the accelerometer input axis approaches the gravity vector. There is then only a small angle between the plane containing the accelerometer input axis (or axes) and the vector to be sensed, this condition being reached whenever the borehole axis approaches horizontal. For the determination of azimuth, the region of poor accuracy is that in which the plane containing the input sensitive axis (or axes) of the acceleration vector sensor and the direction reference vector (either earth rotation or earth magnetic field) approaches parallelism to the plane containing both the earth's gravity vector and the direction reference vector. When using an angular rate sensor and true azimuth is to be computed, this region of poor accuracy exists for a borehole axis that is near true East-true West and near horizontal. When using a magnetic field sensor and magnetic azimuth is to be computed, the region of poor accuracy exists for a borehole axis near magnetic East-magnetic West and near horizontal. In these regions of poor accuracy, small sensor errors will lead to large errors in the desired inclination and/or azimuth.

To overcome these regions of poor accuracy and avoid large errors in such cases, U.S. Pat. Nos. 4,265,028 and 4,197,654 and U.S. patent application Ser. No. 338,261 show that a single vector sensor device can be used to obtain knowledge of three orthogonal components of a reference vector quantity. The method shown in these cited patents and application is that of rotating a sensor about an axis of rotation that has the sensor input sensitive axis canted or inclined relative to a normal to the rotation axis by some angle γ. U.S. Pat. No. 4,265,028 describes the use of a canted accelerometer to measure three orthogonal components of the earth's gravity vector at a fixed (but moveable) location in a well or borehole. U.S. Pat. No. 4,197,654 describes the use of a canted gyroscope (or angular rate sensing device) to measure three orthogonal components of the earth's angular velocity vector. The referenced patent application describes the use of a canted magnetic field sensing device to measure three orthogonal components of the earth's magnetic field vector. Such provision of a cant to the sensor input axis of sensitivity provides a component of the sensed vector along the rotation or borehole axis and this is sufficient to eliminate the geometric regions of poor accuracy which all apparatus having sensing axes only normal to the borehole will have. The present invention discloses apparatus and methods to use this third component of sensed data to be combined with the computation previously described in parent application Ser. No. 351,744 for deriving inclination and azimuth.

SUMMARY OF THE INVENTION

It is a major purpose of this invention to provide method and means to use data from angular rate and acceleration sensors in a mapping or survey tool, one or more of the sensors being canted, to determine the orientation of the inertial angular rate or magnetic field vector sensor (or sensors) with respect to a known earth fixed coodinate set so that correct azimuth determination can be made. It is a second purpose of this invention to provide method and means for azimuth determination in a completely explicit non-ambiguous manner once the sensor data has been resolved to a known earth fixed coordinate set.

The determination of azimuth with respect of either the earth's inertial angular velocity vector (so called true azimuth) or earth's magnetic field vector (so called magnetic azimuth) requires that one first determine at least one (but for complete all azimuths two orthogonal) component of the desired reference vector (angular velocity or magnetic) in a plane parallel to the earth's surface and in a known orientation to the desired unknown azimuth direction. In mapping or survey apparatus of the types cited as previously used in wells or boreholes, the reference direction vector sensors, either inertial angular rate or magnetic, provide outputs proportional to the vector dot product of vectors along their input sensitive axes and the reference vectors. Such outputs of themselves provide no means to know the components of the reference direction vector in a horizontal plane. However, an acceleration sensor (sensors) at a fixed location in the well or borehole provides direct knowledge of the relation of its input axis of sensitivity with respect to the local gravity vector which by definition is normal to the horizontal plane. Since the orientation of the input axis of sensitivity of the reference direction vector sensor, either angular rate or magnetic, is known with respect to the input axis of sensitivity of the acceleration sensor, the output of the acceleration sensor (or sensors) thus may be used to process the output (or outputs) of the reference direction vector sensor (or sensors) to determine one or more components in the horizontal plane.

When the reference direction vector sensor is canted to sense a component along the borehole axis, the third component can be used in computation of azimuth along with the previously stated components resolvable into the horizontal plane. For example, when the acceleration sensor is canted, the component of gravity along the borehole axis may be used in computation with the gravity component in the vertical plane to compute improved accuracy values for the tilt or inclination angle.

Accordingly, it is a major object of the invention to provide borehole survey apparatus wherein angular rate sensor means and acceleration sensor means are suspended and effectively rotated in a borehole, at least one of the sensors may be canted at one angle γ relative to the borehole axis, the angular rate sensor means having amplitude output GA and rotation related phase output GP, and the acceleration sensor means having amplitude output AA and rotation related phase output AP, there also being means supplying a signal vaue Ω_(v) proportional to the local vertical component of the earth's angular rate of rotation, and there being means supplying a value derived from γ, the improvement which comprises

(a) first means for combining AA, AP, GA, GP, said value derived from γ, and Ω_(v) to derive a value ψ for borehole azimuth at the level of the sensor means in the borehole.

In addition, the invention provides means operatively connected with said first means for employing AA modified by the value derived from γ to derive a value θ for borehole tilt from vertical at the level of said sensor means in the borehole.

The basis method of the invention involves the method of borehole mapping or surveying typically using a single angular rate sensor and a single acceleration sensor, both with input axis of sensitivity, one or both sensors being typically canted, the sensors being effectively rotated about the borehole axis, the sensors having outputs.

The outputs of the angular rate sensor and the acceleration sensor are typically employed to derive, from the rate sensor, two or three components respectively in a horizontal plane, one normal to the plane containing the borehole axis and the gravity vector, and the other two or three in that plane, borehole azimuth being derived form the components in a horizontal plane. Three components are formed when the sensors are canted.

These and other objects and advantages of the invention, as well as the details of an illustrative embodiment, will be more fully understood from the following specification and drawings, in which:

DRAWING DESCRIPTION

FIG. 1 is a geometrical depiction of a reference coordinate system established at the start of borehole drilling;

FIG. 1a relates the FIG. 1 co-ordinate system to an instrument level co-ordinate system in a borehole;

FIGS. 2 and 2a show plots of single axis accelerometer and gyroscope outputs vs instrument rotation angle;

FIG. 3 is a geometrical showing of vector relationships in a borehole;

FIG. 4 is a circuit block diagram;

FIG. 5 is a coordinate system diagram;

FIG. 6 is a circuit block diagram;

FIG. 7 shows instrumentation in a borehole (single axis angular rate sensor, and single axis accelerometer);

FIG. 8 shows instrumentation in a borehole (dual axis angular rate sensor, and dual axis accelerometer);

FIG. 9 is an elevation taken in section to show one form of instrumentation employing the invention;

FIG. 10 is an elevation showing use of the FIG. 9 instrumentation in multiple modes, in a borehole;

FIG. 11 is a vertical section showing further details of the FIG. 9 apparatus as used in a borehole; and

FIG. 12 is a circuit block diagram.

DETAILED DESCRIPTION

Referring first to FIG. 9, a carrier such as elongated housing 10 is movable in a borehole indicated at 11, the hole being cased at 11a. Means such as a cable to travel the carrier lengthwise in the hole is indicated at 12. A motor or other manipulatory drive means 13 is carried by and within the carrier, and its rotary output shaft 14 is shown as connected at 15 to an angular rate sensor means 16. The shaft may be extended at 14a, 14b and 14c for connection to first acceleration sensor means 17, second acceleration sensor means 18, and a resolver 19. The accelerometers 17 and 18 can together be considered as means for sensing tilt. These devices have terminals 16a - - - 19a connected via suitable slip rings with circuitry indicated at 29 carried within the carrier (or at the well surface, if desired).

The apparatus operates for example as described in U.S. Pat. No. 3,753,296 and as described above to determine the azimuthal direction of tilt of the borehole at a first location in the borehole. See for example first location indicated at 27 in FIG. 2. Other U.S. patents describing such operation are U.S. Pat. Nos. 4,199,869, 4,192,077 and 4,197,654. During such operation, the motor 13 rotates the sensor 16 and the accelerometers either continuously, or incrementally.

The angular rate sensor 16 may for example take the form of one or more of the following known devices, but is not limited to them:

1. Single degree of freedom rate gyroscope

2. Tuned rotor rate gyroscope

3. Two axis rate gyroscope

4. Nuclear spin rate gyroscope

5. Sonic rate gyroscope

6. Vibrating rate gyroscope

7. Jet stream rate gyroscope

8. Rotating angular accelerometer

9. Integrating angular accelerometer

10. Differential position gyroscopes and platforms

11. Laser gyroscope

12. Combination rate gyroscope and linear accelerometer

Each such device may be characterized as having a "sensitive" axis, (or axes) which is the axis about which rotation occurs to produce an output which is a measure of rate-of-turn, or angular rate ω. That value may have components ω₁, ω₂ and ω₃ in a three axis co-ordinate system. The sensitive axis may be generally normal to the axis 20 of instrument travel in the borehole, or canted at an angle γ.

The acceleration sensor means 17 may for example take the form of one or more of the following known devices; however, the term "acceleration sensor means" is not limited to such devices:

1. one or more single axis accelerometers

2. one or more dual axis accelerometers

3. one or more triple axis accelerometers

Examples of acceleration sensors include the accelerometers disclosed in U.S. Pat. Nos. 3,753,296 and 4,199,869, having the functions disclosed therein. Such sensors may be supported to be orthogonal to the carrier axis. They may be stationary or carouseled, or may be otherwise manipulated, to enhance accuracy and/or gain an added axis or axes of sensitivity. The axis of sensitivity is the axis along which acceleration measurement occurs.

FIG. 11 shows in detail dual input axis rate sensor means and dual output axis accelerometer means, and associated surface apparatus. In FIG. 11, well tubing 110 extends downwardly in a well 111, which may or may not be cased. Extending within the tubing is a well mapping instrument or apparatus 112 for determining the direction of tilt, from vertical, of the well or borehole. Such apparatus may readily be traveled up and down in the well, as by lifting and lowering of a cable 113 attached to the top 114 of the instrument. The upper end of the cable is turned at 115 and spooled at 116, where a suitable metter 117 may record the length of cable extending downwardly in the well, for logging purposes.

The apparatus 112 is shown to include a generally vertically elongated tubular housing or carried 118 of diameter less than that of the tubing bore, so that well fluid in the tubing may readily pass, relatively, the instrument as it is lowered in the tubing. Also, the lower terminal of the housing may be tapered at 119, for assisting downward travel or penetration of the instrument through well liquid in the tubing. The carrier 118 supports first and second angular sensors such as a rate gyroscopes G₁ and G₂, and accelerometers 120 and 121, and drive means 122 to rotate the latter, for travel lengthwise in the well. Bowed springs 170 on the carrier center it in the tubing 110.

The drive means 122 may include an electric motor and speed reducer functioning to rotate a shaft 123 relatively slowly about a common axis 124 which is generally parallel to the length axis of the tubular carrier, i.e. axis 124 is vertical when the instrument is vertical, and axis 124 is tilted at the same angle form vertical as is the instrument when the latter bears sidewardly against the bore of the tubing 110 when such tubing assumes the same tilt angle due to borehole tilt from vertical. Merely as illustrative, for the continuous rotation case, the rate of rotation of shaft 124 may be within the range 0.5 RPM to 5 RPM. The motor and housing may be considered as within the scope of means to support and rotate the gyroscope and accelerometers.

Due to rotation of the shaft 123, and lower extensions 123a, 123b and 123c thereof, the frames 125 and 225 of the gyroscopes and the frames 126 and 226 of the accelerometers are typically all rotated simultaneously about axis 124, within and relative to the sealed housing 118. The signal outputs of the gyroscopes and accelerometers are transmitted via terminals at suitable slip ring structures 125a, 225a, 126a and 226a, and via cables 127, 127a, 128 and 128a, to the processing circuitry at 129 within the instrument, such circuitry for example including that to be described, and multiplexing means if desired. The multiplexed or nonmultiplexed output from such circuitry is transmitted via a lead in cable 113 to a surface recorder, as for example include pens 131-134 of a strip chart recorder 135, whose advancement may be synchronized with the lowering of the instrument in the well. The drivers 131a - - - 134a for recorder pens 131-134 are calibrated to indicate borehole azimuth, degree of tilt and depth, respectively, and another strip chart indicating borehole depth along its length may be employed, if desired. The recorder can be located at the instrument for subsequent retrieval and read-out after the instrument is pulled from the hole.

The angular rate sensor 16 may take the form of gyroscope G₁ or G₂, or their combination, as described in U.S. Pat. No. 4,199,869. Accelerometers 126 and 226 correspond to 17 and 18 in FIG. 9.

Consider now a reference coordinate system is established at the start of the borehole such that X is parallel to the earth surface and North, Y is parallel to the earth surface and East, and Z is perpendicular to the earth surface and Down. The starting point is at latitude λ and FIG. 1 shows the basic geometry.

From this starting reference, the bore axis is defined as rotated through an azimuth angle ψ clockwise about Z, followed by a rotation φ about the new Y axis to obtain a bore axis reference coordinate set in the bore such that z is downward along the bore axis, y is parallel to the earth surface, and x lies perpendicular to y and z. Also, as will be seen, x is in the vertical plane containing the gravity vector and the borehole axis z. See also FIG. 1a.

It may be shown that the direction cosine matrix relating this bore axis reference set to the starting reference set is as shown below:

    ______________________________________                                                       co-ord set at surface                                                          .sup.-- .sub.--X                                                                       .sup.-- .sub.--Y                                                                         Z                                              ______________________________________                                         co-ord. set      x      CψCφ                                                                           SψCφ                                                                           -Sφ                                in Bore          y      -Sψ Cψ  0                                      Hole             z      CψXφ                                                                           SψSφ                                                                           Cφ                                 ______________________________________                                    

In the above C represents Cosine, S represents Sine.

There is no direct way to measure all three direction cosines relating z to the fixed (starting) reference set. However, gyroscopic and acceleration sensing devices can be used to sense quantities from which the required coefficients can be calculated.

The earth rate rotation vector, Ω, in reference coordinates is ##EQU2## where Ω_(H) =ΩCλ, the horizontal compoment,

Ω_(v) =ΩSλ, the vertical component,

and the earth's gravity vector, g, in reference coordinates is ##EQU3##

In the above expression the symbol 1_(n) is a unit vector in the N direction.

The components of Ω and gin the bore axis reference set can be found by forming the dot products 1_(x) ·Ω, 1_(y) ·Ω,1_(z) ·Ω and 1_(x) ·g, 1_(y) ·g, and 1_(z) ·g The results of these operations are: ##EQU4##

Ideally, three accelerometers and three gyro sensing axes could determine all of the required information with no ambiguities or unusual sensitivities other than the classical and well known increased sensitivity to gyro error as latitude increases toward the polar axis.

As shown by the earlier cited patents, to reduce the size and system complexity, sufficient information may be obtained by either a single axis gyro and single axis accelerometer rotated such that their input axes are swept about the x, y plane (normal to the bore axis) or by a two axis gyro and two axis accelerometer having their input axes in the x, y plane. FIG. 7 shows such a single axis gyro G and single axis accelerometer A (see also 16 and 17 in FIG. 9); and FIG. 8 shows such a two axis gyro G₁ and G₂ and two axis accelerometer A₁ and A₂ (see also FIG. 11).

If single axis instruments are used, the plot of their outputs vs rotation angle will appear as (in the absence of sensor errors) shown in FIG. 2.

In this figure, it is evident that the accelerometer output is a sinusoid having its peak output at the point where the input sensitive axis is parallel to the x axis, where x was as previously defined to be in the vertical plane containing the gravity vector and the borehole axis. If the phase angle α, between the accelerometer peak output and the gyro peak output is measured by suitable signal processing, it is then possible to compute Ω_(x), the component of the earth rotation vector in the vertical plane containing the borehole axis and the earth's gravity vector, and Ω_(y), the component of the earth rotation vector in the horizontal plane (normal to the gravity vector). Such components are: ##EQU5## From the previously shown mechanization that Ω_(y) was equivalent to:

    Ω.sub.y =-Ω.sub.H sin ψ                    (11)

it would be possible to compute azimuth as: ##EQU6## using the value of Ω_(y) computed from the gyro output and the phase angle between the gyroscope output and the accelerometer output. This displays the essential usage of the accelerometer output to determine a component of the earth's inertial angular rate vector in a horizontal plane.

The method shown above is suitable except that for azimuths near 90° (East) or 270° (West) the arcsin function provides very poor sensitivity since the rate of change of sin ψ with ψ is very low in these regions. This would lead to large errors in azimuth from small sensor errors. It is, therefore, desirable to find another component of the earth's inertial angular rotation vector in the horizontal plane. FIG. 3 shows a side view of the borehole along the previously defined y axis.

The value of the measured component Ω_(x) is by inspection:

    Ω.sub.x =Ω.sub.v sin φ+Ω.sub.B cos φ(13)

Since Ω_(x) has been determined from the gyro output and the accelerometer to gyro phase angle, and since φ can be determined from the amplitude of the accelerometer signal as: ##EQU7## Then ##EQU8## But as previously defined; ##EQU9## so that ##EQU10## From this it is possible to compute: ##EQU11##

This mechanization also provides a value of ψ and since it is based on a arccos vs the previously cited arcsin function, the region of poor sensitivity is near azimuth of 0° (North) or 180° (South). This again shows the essential use of the accelerometer output to properly resolve the gyro output into the horizontal plane. If one desires, these two functions can be combined into one which has no regions of poor sensitivity. Such a form is: ##EQU12##

FIG. 4 shows a complete block diagram of the described mechanization. As the combination of sensing devices is rotated about its rotation axis in a borehole, both the inertial angular rate sensing and acceleration sensing devices will produce variable output indications proportional to the vector dot product of a unit vector along the respective input axis and the local earth rotation vector and gravity vector respectively. For continuous rotation operation at a fixed location in the borehole these signals will be sinusoidal in nature. For discrete step rotation, the sensor outputs will be just the equivalent of sampling points on the above mentioned sinusoidal signals. Thus, from a knowledge of sample point amplitudes and position along the sinusoid, the character of an equivalent sinusoid in amplitude and phase may be determined. For either continuous rotation or discrete positioning, the quantities that must be determined are the gyro signal amplitude GA (GAMPLITUDE), the accelerometer signal amplitude AA (AAMPLITUDE), and the phase angle between the peak values of these two signals, α. FIG. 4 shows the two sensor signals, after required scaling, and a reference time or angle signal as inputs to the two blocks labeled "Sinusoid Amplitude and Phase." Each of these blocks finds the amplitude of the input sinusoid and the phase angle between the input signal and the reference derived from the rotation drive function. The outputs of the upper block are gyro amplitude and phase, labelled GA (GAMPLITUDE) and GP (GPHASE). The outputs of the lower block are accelerometer amplitude and phase, labelled AA (AMPLITUDE) and AP (APHASE). These amplitude functions are then directly input to subsequent elements and the required phase difference α, is shown, GPHASE minus APHASE.

If a two axis gyro and two axis accelerometer are used, allowance must be made for unknown rotation β about the bore axis. FIG. 5 shows a view looking at the x, y plane from the positive z side.

The sensed accelerometer outputs in terms of the gravity components g_(x) =g sin φ and g_(y) =0 are:

    g.sub.1 =g.sub.x cos β                                (21)

    g.sub.2 =-g.sub.x sin β                               (22)

The sensed angular rates for the two gyro outputs are:

    Ω.sub.1 =Ω.sub.x cos β+Ω.sub.y sin β(23)

    Ω.sub.2 =-Ω.sub.x sin β+Ω.sub.y cos β(24)

From the accelerometer data ##EQU13##

Using the value of β determined above from the accelerometer data and the sensed outputs of the two gyros, two components of the gyro output Ω_(x) and Ω_(y) in a known coordinate set may be computed as:

    Ω.sub.x =Ω.sub.1 cos β-Ω.sub.2 sin β(27)

    Ω.sub.y =Ω.sub.1 sin β+Ω.sub.2 cos β(28)

These values are then identical to the AAMPLITUDE, Ω_(x), and Ω_(y) previously described for the rotated single axis case and the value of azimuth may be determined in the same way. FIG. 6 shows a complete block diagram of circuitry to perform this determination. In FIG. 6, the differences in signal processing for the two angular rate inputs and the two acceleration inputs compared to one rotated sensor of each kind are as shown in FIG. 4. The portion to the left of the dotted line in FIG. 6 would be substituted for the corresponding portion of FIG. 4. Note that since there are two nominally orthogonal signals of each kind, no reference time or angle input is required. Again, the essential use of acceleration sensor outputs to resolve the angular rate sensor data to a known coordinate system is shown.

Although the previous description has used the earth's inertial angular rate vector as the reference direction vector, the earth's magnetic field vector can be used as a reference if magnetic vector sensors replace angular velocity sensors, as in the drawings. All that is necessary is to substitute M_(H) and M_(V) for Ω_(H) and Ω_(v), M_(x) and M_(y) for Ω_(x) and Ω_(y), and M_(B) for Ω_(B). In these formulations the various components of the earth's magnetic field vector are used and the resulting azimuth is the magnetic azimuth measured with respect to the horizontal component of the earth's magnetic field. The same essential dependence on the acceleration sensors, for the resolution of the magnetic sensor outputs into a horizontal plane, is evident in this usage.

Referring now in detail to FIG. 4, the angular rate sensor (gyroscope) amplitude and phase outputs are indicated at GA and GP. These are typically in voltage signal form. Similarly, the accelerometer amplitude and phase outputs are indicated at AA and AP. A synchronizing reference time or angle signal is supplied at 150 to the amplitude and phase detectors 148 and 149 which respond to the gyroscope and accelerometer outputs to produce GA, GP, AP and AA. Means is also provided to supply at 151 a signal corresponding to earth's rotation rate Ω, and to supply at 152 a signal corresponding to the borehole latitude λ. A sin/cos generator 153 operates on signal 152 to produce the output sin λ at 154. The latter and signal 151 are supplied to multiplier 155 whose output Ω sin λ=Ω_(v) appears on lead 156.

In accordance with the invention, (a) first means is provided for combining (or operating upon) AA, AP, GA, GP and Ω_(v) to derive a value ψ for borehole azimuth at the level of the sensors suspended in the borehole. The azimuth signal ψ appears on lead 157 at the right of the circuitry shown. In addition, (b) second means is operatively connected with the referenced first means for employing AA to derive a signal value φ representative of borehole tilt from vertical, at the level of the sensor means in the borehole.

More specifically, such (a) first means include (c) means responsive to GA, GP and AP to derive:

(i) a first component Ω_(x) of the angular rate sensor output, and

(ii) a second component Ω_(y) of the angular rate sensor output. See Ω_(v) on lead 158, and Ω_(y) on lead 159. Such (c) means may typically include:

(d) means responsive to GP and AP to produce a phase angle value or signal α representative of the difference in phase of the GP and AP signals (see for example the subtractor 159 connected with the output sides of 148 and 149, the subtractor output α appearing on lead 160),

(e) means responsive to α to produce signal values sin α and cos α (see for example the sin/cos generator 161 whose input side is connected with lead 160, and whose outputs sin α and cos α appear on leads 162 ad 163),

(f) means responsive to GA and cos α to multiply same and produce the signal value Ω_(x) (see for example the multiplier 164 whose inputs are connected with GA lead 165 and cos α lead 163),

(g) means responsive to GA and sin α to multiply same and produce the signal value Ω_(y) (see for example multiplier 166 whose inputs are connected with the GA input and with sin α lead 112).

The (a) first means also includes (h) means responsive to Ω_(x), AA and Ω_(v) to derive a value Ω_(B) and (j) means responsive to Ω_(y) and Ω_(B) to derive the said value ψ for borehole azimuth. For example, the (h) means may include:

(h₁) an arcsin generator 170 responsive to AA (supplied on lead 171 from detector 149) to generate output at 172,

(h₂) sin/cos generator 173 connected with output 172 to produce output sin φ, on lead 174 and output cos φ on lead 175,

(h₃) multiplier 176 responsive to sin φ on lead 174 and Ω_(v) on lead 156 to produce their product on output lead 177,

(h₄) subtractor 178 connected with leads 177 and 158 to produce the value (Ω_(x) -Ω_(v) sin φ) on lead 179

(h₅) a divider 180 to divide the values on leads 179 and 175 and produce the desired values Ω_(B) on lead 181.

The (i) means referred to above is shown in FIG. 4 to include an arc tangent generator 182 connected with leads 159 and 181 to be responsive to Ω_(y) and Ω_(B) to produce the ψ output proportional to arctan ##EQU14## The tilt output signal φ is produced on lead 184 connected with the output of arcsin generator 170.

FIG. 6 shows similar connections and circuit elements responsive to inputs Ω₁ and Ω₂ from two gyroscopes (or dual axis gyroscope) and inputs ₁ and ₂ from two accelerometers (or from a dual axis accelerometer), to produce the values Ω_(y), Ω_(x) and AA, which are then processed as in FIG. 4. See also FIG. 8.

The above operational devices as at 148, 149, 159, 178, 155, 164, 166, 180, 153, 161, 173, 170 and 182 may be analog or digital devices, or combinations thereof.

Referring to FIG. 10, the determinations of azimuth ψ and tilt φ are carried out at multiple locations in a borehole, as at 27, 27' and 27"; and they may be carried out at each such location during cessation of elevation or lowering by operation of cable 12, or during such elevation or lowering.

When canted sensors are used, the computations are modified to use the components of the reference direction and gravity vectors along the borehole axis. For the case of single axis sensors, FIG. 2 would be modified to incorporate steady outputs for both the gyro and accelerometer. These appear at 210 and 211 in FIG. 2a. Each steady component, in the absence of sensor errors, is proportional to the component along the borehole rotation axis. Also, the amplitudes of the sinusoids would be reduced. Specifically, if the canted angle is designated as gamma, γ, then as shown in FIG. 2a, ##EQU15##

If the values of AAMPLITUDE and GAMPLITUDE obtained from the canted sensors are modified by a value derived from γ, as for example divided by cos γ, then the previously discussed values g_(x) and √Ω_(x) ² +Ω_(y) ² are obtained and computations can proceed as previously described to compute the resolved gyro components Ω_(x) and Ω_(z), the tilt or inclination angle φ, and the azimuth ψ by any of the three indicated methods. See for example in FIG. 4 the optional provision of the cant angle γ signal source 190, the cos γ generator 191, divider 192 to divide GA by cos γ, and the quotient GA' (which is the modified GA) on lead 193. Associated switches are shown. Similarly, in FIG. 4, AA is divided at 194 by cos γ to produce modified AA' on lead 195 to produce AA' (the modified AA).

When the values of AAVERAGE and GAVERAGE are divided by sin γ, then g_(z) and Ω_(z) are obtained. Since as previously shown:

    g.sub.x =g sin φ                                       (33)

    g.sub.z =-g cos φ                                      (34)

a value of φ may be computed either as: ##EQU16##

Either of these methods or values is free of the reduced sensitivity for near horizontal boreholes of the previously (FIG. 4) shown: ##EQU17##

Of the two new forms shown, the Arctan form of Equation (36) has the additional benefit that it is not in error due to either a scale factor error, or an error in the knowledge of gravity. Also, as previously shown:

    Ω.sub.z =Ω.sub.H cos ψ sin φ-Ω.sub.v cos φ(8)

Returning to FIG. 3, it is again possible to compute the shown vector Ω_(B) from Ω_(z) rather than from Ω_(x). By inspection

    Ω.sub.z =Ω'B sin φ-Ω.sub.v cos φ (38)

Solving for Ω'_(B) : ##EQU18## But as previously defined,

    Ω.sub.z =Ω.sub.H cos ψ sin φ-Ω.sub.v cos φ(40)

So that:

    Ω'.sub.B =Ω.sub.H cos ψ                    (41)

From this it would be possible to compute azimuth as: ##EQU19## or as previously shown, using Ω_(y) and Ω'_(B), ##EQU20## FIG. 12 shows a block diagram of electrical apparatus for performing these computations.

Referring to detail to FIG. 12, elements thereof also found in FIG. 4 carry the same identifying numerals. Also shown is means supplying a value derived from γ, the cant angle of the gyroscope and accelerometer relative to the borehole axis. Such means is shown for example to include a voltage at 210 corresponding to cant angle γ, and a sine/cosine generator 211 to produce sine γ and cosine γ outputs at 212 and 213. These are, of course, trigonometric values drived from γ.

In accordance with the invention, (a) first means is provided for combining AA, AP, GA, GP, a value or values derived from γ, (i.e. sin γ and cos γ, for example) and Ω_(v) to derive a value ψ for borehole azimuth at the level of the sensor means in the borehole. In addition, (b) second means is operatively connected with the first means for employing AA modified by a value derived from γ to derive a value φ for borehole tilt from vertical at the sensor level in the hole.

More specifically, the first means include (c) means such as divider 215 responsive to GA and cos γ to divide GA by cos γ and produce a value GA', indicated at lead 216. Also, the first means (c) typically is responsive to GP and AP to derive a primary component Ω_(y) of the angular rate sensor output. See for example sin α generator 161 and multiplier 217 for multiplying sin α and GA' to produce Ω_(y), to be used in the derivation of ψ.

The first means is also shown as including means responsive to GAV output 217 from device 148, divided by sin γ at 218 to produce a value Ω_(z) on lead 217. Multiplier 220 then multiplies Ω_(z) by the value Ω_(v) cos φ to produce an output on lead 221, the latter then being divided by sin φ at 222 to produce Ω'_(B). Generator 223 then responds to Ω'_(B) and Ω_(y) to produce ψ by computing arctan Ω_(y) /Ω'_(B), as shown in FIG. 12.

In the above, φ is the borehole tilt angle, and is computed by the devices shown. Thus, output AA at 171 is divided at 230 by cos α to produce g_(x) ; AAV output from device 149 is divided by sin γ at 231 to produce output g_(z), and arctangent generator 224 receives g_(x) and g_(y) to compute arctan -g_(x) /g_(z), which produces φ. Sin/cosine generator 232 receives φ to produce sin φ at 233, and cos φ at 234, for use as described above. Multiplier 235 multiplies cos φ and Ω_(v) to produce Ω_(v) cos φ on line 237.

It is of course possible to have more than one axis of sensing of each kind, canted. If, in the previously described approach using a two axis gyro and a two axis accelerometer, both sensor axes of each kind are canted, then two independent estimates of the component of the sensed reference vector along the borehole axis are obtained.

Although the previously described computations used the example of a canted accelerometer and a canted angular rate sensor, it is clear that either sensor could be canted alone or that a magnetic field sensor could be substituted for the angular rate sensor if magnetic azimuth were desired. 

We claim:
 1. In mapping or survey apparatus using a single inertial angular rate sensor and a single acceleration sensor, one or both with input axes of sensitivity nominally canted relative to a well or borehole axis, rotated about the borehole axis and both sensors having outputs, the combination comprising(a) means responsive to the accleration sensor output together with the inertial angular rate sensor output to compute from the rate sensor output the components thereof in the horizontal plane and the vertical plane, and (b) means to derive azimuth as the Arcsin of the component of sensed inertial angular rate in the horizontal plane normal to the plane containing the borehole axis and the gravity vector divided by the horizontal plane component of the earth's angular velocity vector.
 2. The apparatus of claim 1 including(c) means responsive to the acceleration sensor output to compute the component of the earth's rotation rate in the horizontal plane at its intersection with the vertical plane containing the gravity vector and the borehole axis from the sensed angular rate, and (d) means to derive azimuth as the Arcos of the component computed in (c) divided by the horizontal plane component of the earth's angular velocity vector.
 3. The apparatus of claim 2 including(e) means to compute azimuth as the Arctan of the argument computed in (b) divided by the argument computed in (c).
 4. In mapping or survey apparatus using two axis inertial angular rate sensor means and two-axis acceleration sensor means both with their two input axes of sensitivity nominally normal to a well or borehole axis, both sensors having outputs,(a) means responsive to acceleration sensor outputs together with the inertial angular rate sensor outputs to compute from the rate sensor outputs the components thereof in the horizontal and vertical plane, and (b) means to derive azimuth as the Arcsin of the component of sensed inertial angular rate in the horizontal plane normal to the plane containing the borehole axis and the gravity vector divided by the horizontal plane component of the earth's angular velocity vector.
 5. The apparatus of claim 4 including(c) means responsive to the acceleration sensor output to compute the component of the earth's rotation rate in the horizontal plane at its intersection with the vertical plane containing the gravity vector and the borehole axis from the sensed angular rate, and (d) means to compute azimuth as the Arccos of the component computed in (c) divided by the horizontal plane component at the earth's angular velocity vector.
 6. The apparatus of claim 5 including(e) means to compute azimuth as the Arctan of the argument computed in (b) divided by the argument computed in (c).
 7. The combination of claim 1 where the inertial angular rate sensor is replaced by a magnetic field vector sensor, the earth's angular velocity vector is replaced by the earth's magnetic field vector, and magnetic azimuth is computed.
 8. The combination of claim 1 wherein the inertial angular rate sensors are replaced by magnetic field vector sensors, the earth's angular velocity vector is replaced by the earth's magnetic field vector, and magnetic azimuth is computed.
 9. The apparatus of claim 2 in which one of the sensors has its input axis of sensitivity canted so as to provide sensing of a component along the borehole axis, and employing a canted accelerometer output component along the borehole axis to increase the accuracy of tilt or inclination angle for near horizontal boreholes.
 10. The apparatus of claim 9 includes means to compute tilt or inclination as the Arccos of the gravity component along the borehole axis, divided by the known magnitude of the earth's gravity vector.
 11. The apparatus of claim 9 including means to compute tilt or azimuth as the Arctan of the gravity component normal to the borehole axis in the vertical plane divided by the component along the borehole axis.
 12. The apparatus of claim 9 including means responsive to the acceleration sensor output together with the canted inertial angular rate sensor output component along the borehole axis to compute a second estimate of the component of the earth's rotation rate in the horizontal plane at its intersection with the vertical plane containing the gravity vector and the borehole axis.
 13. The apparatus of claim 9 including means to compute azimuth as the Arccos of the component computed in (d) divided by the horizontal plane component of the earth's angular velocity vector.
 14. The apparatus of claim 9 including means to compute azimuth as the Arctan of the argument computed in (b) of claim 1 above divided by the argument computed in (d) above.
 15. In borehole survey apparatus wherein angular rate sensor means and acceleration sensor means are suspended and effectively rotated in a borehole, at least one of said sensors canted at an angle γ relative to the borehole axis, the angular rate sensor means having amplitude output GA, and rotation related phase output GP, and the acceleration sensor means having amplitude output AA and rotation related phase output AP, there also being means supplying a signal value Ω_(v) proportional to earth's angular rate of rotation, and means supplying a value derived from γ, the improvement which comprises(a) first means for combining AA, AP, GA, GP, said value derived from γ, and Ω_(v), to derive a value ψ for borehole azimuth at the level of said sensor means in the borehole.
 16. The apparatus of claim 15 including(b) second means operatively connected with said first means for employing AA modified by said value derived from γ to derive a value φ for borehole tilt from vertical at the level of said sensor means in the borehole.
 17. The apparatus of claim 15 wherein said first means includes (c) means responsive to GA modified by said value derived from γ, GP and AP to derive(i) a first component Ω_(x) of the angular rate sensor output, and (ii) a second component Ω_(y) of the angular rate sensor output.
 18. The apparatus of claim 17 wherein said (c) means to derive Ω_(x) and Ω_(y) includes (d) means responsive to GP and AP to produce a phase angle value α representative of the difference in phase of said GP and AP outputs, (e) means responsive to α to produce sin α and cos α values, (f) means to multiply modified GA and said sin α value to produce Ω_(y), and (g) means to multiply modified GA and said cos α value to produce Ω_(x).
 19. The apparatus of claim 17 wherein said first means includes (h) means responsive to Ω_(x), AA modified by said value derived from γ and Ω_(v) to derive a value Ω_(B) ', and (j) means responsive to Ω_(y) and Ω'_(B) to derive said value ψ for borehole azimuth.
 20. The apparatus of claim 19 wherein said (h) means includes:(h₁) an arc sin generator responsive to modified AA to generate an output, (h₂) sin and cos generator means responsive to said output of the arc sin generator to generate an output sin φ and an output cos φ, (h₃) multiplier means responsive to sin φ and Ω_(v) to produce a product thereof, (h₄) substractor means responsive to said product and Ω_(x) to obtain a difference value, (h₅) divider means to divide said difference value by said output cos φ to obtain said value Ω'_(B).
 21. The apparatus of either one of claims 19 and 20 wherein said (i) means includes an arc tangent generator responsive to Ω_(y) and Ω_(B) to produce an output proportional to arc tan -Ω_(y) /Ω'_(B) which is representative of azimuth ψ.
 22. The apparatus of claim 20 wherein said elements (h₁)-(h₅) are operatively interconnected.
 23. In well bore survey apparatus wherein angular rate sensor means and accelerometer means are located in a borehole and at a cant angle γ relative to the borehole axis, the angular rate sensor means having amplitude output GA and phase output GP, and the accelerometer means having amplitude output AA and phase output AP, there being means providing a value Ω_(v) proportional to earth's angular velocity vector, and means supplying the combination comprising values derived from γ,(a) means operatively connection to said sensors to be responsive to GA, GP and AP and a value derived from γ to derive a first component Ω_(x) of the angular rate sensor output, (b) means operatively connected to said sensors to be responsive to GA, GP and AP and a value derived from γ to derive a second component Ω_(y) of the angular rate sensor output, (c) means operatively connected to said (a) means to be responsive to Ω_(x), AA and Ω_(v) to derive a value Ω'_(B), and (d) means operatively connected to said (b) and (c) means to derive ψ from Ω_(y) and Ω'_(B) wherein ψ is an azimuth value indicative of the azimuth angle of the borehole relative to the true North at the location of said sensor means.
 24. The combination of claim 23 including(a) means responsive to AA to derive a value φ for borehole tilt at the location of said sensor means in the borehole.
 25. The apparatus of either one of claims 15 and 23 including means suspending said rate sensor means and accelerometer sensor means in the borehole at an elevation at which said derivation of ψ is carried out.
 26. In borehole survey apparatus wherein magnetic sensor means and acceleration sensor means are suspended and effectively rotated in a borehole and at a cant angle γ relative to the borehole axis, the magnetic sensor means having amplitude output GA and rotation related phase output GP, and the acceleration sensor means having amplitude output AA and rotation related phase output AP, there also being means supplying a signal value Ω_(v) proportional to earth's angular rate of rotation, and means supplying a value derived from γ the improvement which comprises(a) first means for combining AA, AP, GA, GP said value derived from γ and Ω_(v) to derive a value ψ for borehole azimuth at the level of said sensor means in the borehole.
 27. The apparatus of claim 26 including(b) second means operatively connected with said first means for employing AA modified by said value derived from γ to derive a value φ for borehole tilt from vertical at the level of said sensor means in the borehole.
 28. The apparatus of claim 26 wherein said first means includes (c) means responsive to GA modified by said value derived from γ, GP and AP to derive(i) a first component Ω_(x) of the magnetic sensor output, and (ii) a second component Ω_(y) of the magnetic sensor output.
 29. The apparatus of claim 28 wherein said (c) means to derive Ω_(x) and Ω_(y) includes (d) means responsive to GP and AP to produce a phase angle value α representative of the difference in phase of said modified GP and AP outputs, (e) means responsive α to to produce sin α and cos α values, (f) means to multiply modified GA and said sin α value to produce Ω_(y), and (g) means to multiply modified GA and said cos α value to produce Ω_(x).
 30. The apparatus of claim 28 wherein said first means includes (h) means responsive to Ω_(x), modified AA and Ω_(x) to derive a value Ω'_(B) and (j) means responsive to Ω_(y) and Ω'_(B) to derive said value ψ for borehole azimuth.
 31. The apparatus of claim 30 wherein said (h) means includes:(h₁) an arc sin generator responsive to modified AA to generate an output, (h₂) sin and cos generator means responsive to said output of the arc sin generator to generate an output sin φ and an output cos φ, (h₃) multiplier means responsive to sin φ and Ω_(v) to produce a product thereof, (h₄) substractor means responsive to said product and Ω_(x) to obtain a difference value, (h₅) divider means to divide said difference value by said output cos φ to obtain said value Ω'_(B).
 32. The apparatus of one of claims 15, 23 and 26 wherein said value derived from γ is a value for cos γ and said first means (a) includes means for dividing GA by said cos γ value to derive a modified value of GA, and for dividing AA by said cos γ value to derive a modified vlaue of AA.
 33. The apparatus of claim 16 wherein said angular rate sensor has a steady output component GAV and said acceleration sensor means has a steady output component AAV, and said first means receives said steady outputs for combination with AA, AP, GA, GP, values derived from γ, and Ω_(v) to derive said value ψ.
 34. The apparatus of claim 33 wherein said values derived from γ correspond to cos γ and sin γ, and said first means includes (c) means responsive to GA divided by cos γ to produce a value GA', GP and AP to derive(i) a primary component Ω_(y) of the angular rate sensor output.
 35. The apparatus of claim 34 wherein said (c) means to derive Ω_(y) includes (d) means responsive to GP and AP to produce a phase angle value α representative of the difference in phase of said GP and AP outputs, (e) means responsive to α to produce sin α and cos α values, and (f) means to multiply GA' and said sin α value to produce Ω_(y).
 36. The apparatus of claim 35 wherein said first means includes (g) means responsive to GAV divided by sin γ to produce a value Ω_(z), and including (h) means responsive to said values AA and AAV, and said cos α and sin α values and _(v) to produce a value Ω_(y) cos φ for addition to Ω_(z) and subsequent division by sin φ to produce a value Ω'_(B).
 37. The apparatus of claim 36 including (i) means responsive to Ω_(y) and Ω'_(B) to produce said value ψ.
 38. The apparatus of claim 36 wherein φ is the borehole tilt angle, and said (h) means includes (j) means responsive to AA, AAV, cos α and sin α to generate φ.
 39. The apparatus of claim 15 wherein the angular rate sensor is canted at said angle γ.
 40. The apparatus of claim 15 wherein the acceleration sensor means is canted at said angle γ.
 41. Apparatus as defined in claim 9 wherein the inertial angular rate sensor is replaced by a magnetic field sensor, the earth's angular velocity vector is replaced by the earth's magnetic field vector and magnetic azimuth is computed. 